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Observation of high angular momentum excitons in cuprous oxide 
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The recent observation of dipole-allowed P-excitons up to principal quantum numbers of n = 25 
in cuprous oxide has given insight into exciton states with unprecedented spectral resolution. While 
so far the exciton description as a hydrogen-like complex has been sufficient for cubic crystals, 
we demonstrate here distinct deviations: The breaking of rotational symmetry leads to mixing 
of high angular momentum F- and P-excitons with the P-excitons so that they can be observed 
in absorption. The P-excitons show a three-fold splitting that depends systematically on n, in 
agreement with theoretical considerations. From detailed comparison of experiment and theory we 
determine the cubic anisotropy parameter of the CU 2 O valence band. 
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Introduction. Excitonic effects are decisive for the op¬ 
tical properties of semiconductors and insulators [1]. Not 
only leads the Coulomb interaction between an electron 
and a hole to a series of bound states, the excitons, with 
energies below the band gap, but also above the gap the 
Coulomb effects lead to a massive redistribution of oscil¬ 
lator strength towards the low-energy states compared to 
a free particle description. Due to this importance it has 
been a major goal to develop a detailed understanding 
of excitons on a quantitative level [1]. The description 
of the bound exciton states by the hydrogenic model has 
turned out to be extremely successful in this respect, in 
particular, for bulk semiconductors of cubic symmetry. 

For excitons with wavefunction extensions much larger 
than the crystal unit cell (the Mott-Wannier excitons) 
the hydrogen formula for their binding energy, TZjn^ with 
the Rydberg energy 7^ in a state of principal quantum 
number n, can be simply adapted to the solid state case 
by (i) changing the reduced mass of electron and pro¬ 
ton m to that of electron and hole m*, and (ii) screen¬ 
ing the carrier interaction by the dielectric constant e: 
TZ* = TZm* ! {e^m). The influence of the many-body crys¬ 
tal environment is thus comprised in material properties 
that for cubic semiconductors are, as a rule, isotropic 
such as the scalar dielectric constant e, leading to a for¬ 
mula for excitonic energies that is identical to the one in a 
system with rotational symmetry. The material environ¬ 
ment typically causes a reduction of the atomic Rydberg 
energy by 2 — 3 orders of magnitude into the meV range. 

For the hydrogen problem the spatial symmetry is de¬ 
termined by the continuous rotation group SO (3), where 
the square of the orbital momentum = l{l -\- 1 ) 1 t? 
and its z-component = mh, the magnetic quantum 
number, are constants of motion, making the problem 
integrable. Due to this symmetry the energy levels are 
degenerate with respect to m. For the 1/r-dependence 
of the Coulomb potential, the Lenz-Runge vector is also 
conserved as specific consequence of the underlying SO (4) 
symmetry, causing the energy level degeneracy in 1 . 

The latter degeneracy is lifted for hydrogen-like sys¬ 
tems such as Rydberg atoms [2], where one electron is 


excited into a shell with n ^ 1. Here the screening of the 
nuclear Coulomb potential by inner shell electrons causes 
a deviation from the l/r behavior, breaking the SO(4) 
symmetry. Phenomenologically this can be described by 
the quantum defect model, in which the binding energy 
formula is modified to 7?./(n — 5iY. The quantum defect, 
5i, depends on 1, as the screening varies with the angular 
momentum of the outmost electron state due to its differ¬ 
ent penetration into the electron core. The m-degeneracy 
is still maintained, though, for rotational symmetry. 

In crystals this rotational symmetry is broken down 
to discrete groups. For bulk cubic crystals these are 
the groups Oh or T 4 (depending on presence or ab¬ 
sence of inversion symmetry). Since both groups rep¬ 
resent still quite high symmetry, e.g., compared to low¬ 
dimensional semiconductors, the deviations from the ro¬ 
tational symmetry are usually captured solely by the 
lattice-periodic Bloch functions of electron and hole. The 
envelope wave function of Mott-Wannier excitons, how¬ 
ever, are typically considered as hydrogen-like involving 
the spherical harmonics as angular momentum eigenfunc¬ 
tions, even though, strictly speaking, angular momentum 
is no longer conserved, and l,m are not good quantum 
numbers. So far, no indications for a failure of this de¬ 
scription such as observation of a splitting of levels of 
particular I and/or mixing of levels with different I have 
been reported. 

Due to the small Rydberg energy in prototype semi¬ 
conductors of highest quality like GaAs {TZ* = 4.2 meV), 
consequences of the reduced crystal symmetry are hard to 
resolve in optical spectra of excitons. One might seek for 
deviations from the exciton hydrogen model in materials 
with larger Rydberg energy such as oxides or nitrides, for 
which possible exciton level splittings, however, may be 
blurred by crystal inhomogeneities. A somewhat unique 
position in terms of crystal quality is held by cuprous 
oxide (CU 2 O) with a Rydberg energy of about 90 meV. 
The high quality of CU 2 O natural crystals allowed the 
first experimental demonstration of excitons [3, 4], and it 
is evidenced also by remarkably narrow absorption lines 
as highlighted by the observation of the paraexciton in 
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magnetic field with a record linewidth below 100 neV [5]. 

Here we report a high resolution absorption study of 
the yellow exciton series in cuprous oxide. Besides the 
dominant P-excitons we resolve exciton triplets from 
n = 4 to 10 with linewidths in the ^eV-range. These 
lines are ascribed to P-excitons (more precisely, to mixed 
P- and P-excitons). Their emergence in the spectra and 
their splitting which depends systematically on the prin¬ 
ciple quantum number n are consequences of the reduc¬ 
tion from full rotational to discrete Oh symmetry. These 
findings are in good agreement with calculations in which 
the band structure details are taken into account. 

Symmetry analysis. The yellow exciton series is asso¬ 
ciated with transitions between the highest valence and 
lowest conduction bands in CU 2 O. In this material the 
topmost valence band corresponds to the irreducible rep¬ 
resentation Vh = of the Oh point group, and the 
bottom conduction band corresponds to Pe = T)!", see 
Ref. [2] for details. The exciton wavefunction transforms 
according to the product = Pg x Vh x P^ where P^. 
is the representation describing the wavefunction of the 
electron-hole relative motion [ 8 ]. 

To draw the analogy with hydrogen, it is convenient 
to follow the conventional description and categorize the 
excitonic states by the symmetry of the relative motion 
envelope P^.. For P-excitons and P-excitons P^ is irre¬ 
ducible, P 5 = T)^ and Pp = FJ, respectively, while for 
the states with higher angular momentum, e.g. the F 
{I = 3) and H {I = 5) excitons the representations are 
reducible, Pp = F^ + FJ -|- F^ and Pp = F;^ -f 2 r 4 -|- F^, 
demonstrating that I is no longer a good quantum num¬ 
ber for cubic symmetry. Still, we will maintain the exci¬ 
ton classification according to orbital angular momentum 
for simplicity. 

To be optically active, the exciton state representa¬ 
tion Pa; has to contain the three-dimensional = FJ 
representation, corresponding to the components of the 
electric dipole operator. The product of the electron and 
hole Bloch functions can be decomposed as Pe x Vh = 
F^ +r^. Taking into account the envelope function, the 
resulting P-excitons are optically active in one-photon 
transitions because F^ is contained in the product rep¬ 
resentation: F4 G Pp X F^, while this is not the case for 
the S- and P-excitons, so that they are dark [ 6 ]. The ori¬ 
gin of this behavior is the even parity of conduction and 
valence bands in CU 2 O so that dipole transitions have 
to involve odd envelopes. However, accounting for the 
cubic crystal symmetry makes the F- and P-excitons 
dipole allowed in addition to the P-excitons: One can 
readily check that the product Pp x F^ contains 3 F 4 
and Pp X FJ contains one more F^. Hence P-exciton 
states can give rise to four spectral lines. However, one of 
those arising from the FJ product of the Bloch functions 
is weak [17]. Analogously, the product Pp x (FJ -|- F^) 
contains bFj, hence P-exciton can give rise to five lines. 

Experimentals. To test these predictions, high- 
resolution absorption spectra were recorded by detecting 
the emission of a frequency-stabilized laser with 1 neV 
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Figure 1: Bottom: absorption spectrum of the CU 2 O yel¬ 
low exciton series in the energy range of states with principal 
quantum numbers from n = 4 to n = 9. The top panels 
show close-ups of the high energy flanks of the P-excitons 
with n = 4, 5, 6, and 7, respectively. T = 1.2 K. 


linewidth (corresponding to about 250 kHz band width) 
after transmission through a 30 fam. thick CU 2 O crystal 
slab. The wavelength of the laser emission was scanned 
in the range of interest from 570 to 580 nm, for details 
see Ref. [6]. The sample was held in liquid Helium at 
a temperature of 1.2 K and was strain-free mounted in 
a holder that allowed also application of an electric field 
along the optical axis. 

The bottom panel of Fig. 1 shows an absorption spec¬ 
trum of the yellow exciton series in the energy range cor¬ 
responding to principal quantum numbers from n = 4 
to n = 9. The spectrum is dominated by strong ab¬ 
sorption features of the P-excitons, discussed in detail in 
Ref. [6] . However, on the high energy flank of these lines 
weak additional features appear. We highlight that these 
features can be observed starting from n = 4 only. The 
top panels show close-ups of the high energy flanks of the 
n = 4, 5,6, and 7 P-excitons. Each group of features con¬ 
sists of a triplet of lines with the splitting between them 
decreasing systematically with increasing n. The width 
of each line is in the /ieV-range, also decreasing with n 
which corresponds to lifetimes in the nanosecond range. 
Starting from n = 6 another feature appears on the high 
energy side of the triplet, approaching the triplet with 
increasing principal quantum number. 

For n > 7 the triplet is too close to the P-excitons to be 
resolved in absolute transmission. Therefore we have ap¬ 
plied modulation spectroscopy by recording the differen¬ 
tial absorption with and without an electric field applied. 
Only a small field of 15V/cm is applied to avoid notable 
modifications both of absolute exciton energies as well 
as splittings between them. The bottom panel in Fig. 2 
shows such a modulation spectrum, in which triplets can 
be seen up to higher principal quantum numbers. For 
comparison with Fig. 1, the energy range around the 
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Figure 2: Modulation spectra recorded as difference of trans¬ 
mission spectra with and without a small electric field applied. 
The electric field strength is 15 V/cm. The bottom panel 
shows the spectral range of excitons from n = 7 up to 10; 
the top panels give close-ups of the corresponding F-exciton 
features. 


n = 7 exciton is shown, where one recognizes the triplet 
at the same energies as before. For n = 8, also the triplet 
can be seen, where however the two low energy lines have 
almost merged, while the high energy line is still well sep¬ 
arated from them. For higher n only a doublet of lines 
can be seen, as the splitting between the low energy lines 
is comparable or smaller than their line widths. 

The angular momentum of exciton states belonging to 
a particular n is limited by n — 1, hence, in the hydro¬ 
gen model the F-states appear only for n ^ 4. Since 
the observed triplet indeed emerges in Fig. 1 starting 
from n = 4 only, we therefore assign it to F’-excitons. 
We note that based on the symmetry considerations pre¬ 
sented above, F-excitons belonging to a particular n are 
expected to appear in optical spectra as a quadruplet 
(with one line being weak). The numerical calculations 
presented below confirm this and demonstrate that one of 
the lines is indeed substantially weaker as compared with 
three others. Moreover, its energy is very close to that 
of one strong line. For additional confirmation, experi¬ 
ments in moderate magnetic fields (not shown here) were 
carried out that demonstrate a splitting into 6 lines, in 
agreement with our symmetry analysis which predicts a 
two-fold splitting of each line in the triplet. Furthermore, 
the additional feature emerging from n = 6 onwards can 
be attributed to 7F-excitons. Their expected splitting is, 
however, too small to be resolved. 

To the best of our knowledge, such high angular mo¬ 
mentum exciton states including their fine structure split¬ 
ting have not been resolved so far. We emphasize here 
that the description as a hydrogen-like complex can ex¬ 
plain neither the optical activity of F- and iJ-excitons, 
nor their splitting. It is therefore a unique signature of 
the breaking of the rotational symmetry in the cubic crys¬ 
tal and a consequence of its discrete symmetry. For a 


detailed understanding we have developed a microscopic 
description of the F-exciton fine structure. 

Microscopic theory. In CU 2 O the excitonic Rydberg 
TZ* « 90 meV and the splitting between the Fy and 
F^ valence bands A « 130 meV have the same order of 
magnitude. Therefore, a consistent theory of excitonic 
states has to be based on treating the complex valence 
band structure and the Coulomb interaction on the same 
level. Correspondingly, we follow the approach of Ref. [2] 
and present the exciton Hamiltonian (for zero center-of- 
mass wavevector) in the form 

^ ^ ¥~r~^ 

Here p is the momentum of the relative electron-hole mo¬ 
tion, I is the angular momentum one operator acting in 
the basis of the orbital hole Bloch functions F^, and 
is the hole spin operator {sh = 1/2). The energies are 
measured in units of the “bulk” excitonic Rydberg TZ* = 
e^mo/( 2 h^£^ 7 (), the distances are measured in units of 
the corresponding Bohr radius, a* = h^£ 7 (/(e^mo), e 
is the static dielectric constant, A = A/TZ* is the di¬ 
mensionless splitting between the F^ and F^ bands, 
7i = 7i + rnolnie, /r = (673 -I- 472 )/( 57 (), We is the 
conduction electron mass, and the 7 ^ (z = 1 , 2 ,3) are the 
Luttinger parameters. In contrast to Ref. [2] we include 
into the Hamiltonian ( 8 ) the contribution of the cubic 
symmetry responsible for the valence band warping [ 6 ] 

«= = 4 f E X ■'‘“’I™ + X /<^)]»>) , 

\fe=±4 / 

( 2 ) 

with 6 = (73 — 72 )/ 7 i, This extension is crucial, as our 
calculations show, to describe the fine structure of F- 
exciton states absent otherwise. The central-cell correc¬ 
tion to the Coulomb potential as well as the short-range 
electron-hole exchange interaction are disregarded since 
for P- and F-excitons the wave function of the relative 
motion vanishes for coinciding electron and hole coordi¬ 
nates. In Eqs. ( 8 ) and (15) we use and for the 
second-rank irreducible components of the tensors PiPj 
and lilj, where i,j = x,y,z and pi, li are the Cartesian 
components of p and I, respectively. We note that the 
quartic terms in the dispersion, like pt +Py + Pyj which 
are allowed in Oh, result in a F-F-mixing and make the 
F-states active but do not cause their splitting. 

The Hamiltonian ( 8 ) without cubic contribution "He 
{S = 0) has full rotational symmetry and already provides 
a quite accurate description of the F-exciton state ener¬ 
gies [2]. Therefore, it is instructive to disregard Tic and 
determine the spectrum and wavefunctions of F- and F- 
excitons. Moreover, the suppressed short-range electron- 
hole exchange interaction makes it possible to disregard 
the electron spin and characterize the excitons by the 
hole total momentum F and its z-component F^, where 
F = J-l-X is the sum of the momentum of the hole Bloch 
functions J = Sh + I and the envelope function orbital 
momentum L. 
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Figure 3: Comparison of theoretical calculations and exper¬ 
imental data for P-exciton binding energy (a), splitting be¬ 
tween F- and P-excitons (b), and splittings of P-triplet as 
defined in Fig. 1 (c) vs principal quantum number. Dots give 
experimental data. Theory is shown by lines connecting cal¬ 
culated data for discrete values of n. The parameters of the 
calculation are: TV = 87 meV, /r = 0.47, A = 134 meV [2], 
and 5 = —0.1. 


The P-envelopes of excitons correspond to F = 1/2 
and F = 3/2 states [2], while F-excitons correspond in 
this approximation to F = 5/2 and 7/2, respectively. 
Making use of its symmetry one can express the wave- 
function of a state with given F via basic functions with 
total momentum F, |FJF), and radial functions ghjiv), 
= '^Lj 9 lj{'i’)\LJF) (see Refs. [2, 4] for details), 
where the functions gLjij") are determined numerically 
following the method described in Refs. [4, 12], see Sup¬ 
plement [13] for details. 

With inclusion of the cubic symmetry, F and are 
no longer good quantum numbers. In particular, the 
F = 5/2 state gives rise to F/" and Fg states while 
F = 7/2 gives rise to Fg , F/ and Fg states. Note, that 
Fg and F^ are optically active (F 4 G Fg x Fg, Fg x Fg), 
while the Fy states are dark. States of the same symme¬ 
try have to be treated as close-to-degenerate, see Refs. [ 6 ] 
and [13] for details. We have calculated the energy spec¬ 
trum of the P- and F-excitons for different values of the 
cubic anisotropy parameter S keeping all other parame¬ 
ters known from literature fixed. Our calculations show 
that the reasonably small value of S = —0.1 gives good 
accord with experiment [13]. In this case two optically 
active states out of four have very close energies. More¬ 
over, one of those states, namely, Fg originating from 
F = 7/2 has small oscillator strength[13]. Smaller/larger 
values of S result in too small/large predictions for the 
F-shell splitting, so that the estimate of 6 from the ex¬ 
perimental data is quite accurate. 

Discussion. Figures 3(a) and (b) compare the results 
of calculations of the F-exciton binding energy and the 
splitting between F- and P- excitons as function of prin¬ 
cipal quantum number n (the solid lines connect cal¬ 


culated values for discrete n) with the measured data 
(dots), and serve as reference for the accuracy of the de¬ 
scription of the exciton states by our model. For the F- 
excitons we have calculated the center-of-gravity of the 
F-exciton lines. The comparison shows excellent agree¬ 
ment between theory and experiment, providing confi¬ 
dence that also the geW splittings between the F-excitons 
can be assessed by our model. 

These splittings are shown in Fig. 3(c), where A 13 is 
the energy separation between the outermost levels and 
A 23 is the energy separation between the middle and 
higher energy level, see Fig. 1. Note, that the calcu¬ 
lated energy of the fourth (weakly) active state coin¬ 
cides within the accuracy of several percent with the 
energy of middle state in the triplet. The points with 
4 ^ n ^ 7 were measured in absolute transmission (see 
Fig. 1), while points with n > 8 were taken from mod¬ 
ulation spectroscopy. Both magnitudes of the splittings 
and their dependence on n are in good agreement with 
our calculations. The splitting between the states de¬ 
creases strongly with n, as higher-n states are more ex¬ 
tended in space so that their wave function averages over 
more crystal unit cells and become less sensitive to crystal 
symmetry deviations from the full rotation group. More 
quantitatively, the decrease is due to the fact that the 
states are intermixed by quadratic combinations of the 
momentum operator pipj. The corresponding interaction 
matrix elements scale as the inverse square of the state 
radii. In the hydrogen model (r~^} oc n~^(l -|- 1)“^, indi¬ 
cating that the splittings for higher angular momentum 
states decrease with n and 1. This explains the missing 
splitting of the F-excitons in absorption. 

The agreement of the developed model with the ex¬ 
perimental data allows us to accurately evaluate the va¬ 
lence band parameters in CU 2 O. Taking the static di¬ 
electric constant e = 7.5 [14], the electron effective mass 
rUe = 0.99mo [15], and using our values for TZ*, p and 
6 (see text above and caption to Fig. 3), we obtain the 
following values of the Luttinger parameters: 71 = 1.79, 

72 = 0.82, and 73 = 0.54. Interestingly, in this material 

73 < 72 , which is also indicated by microscopic calcula¬ 
tions [16], see Supplement. We emphasize that the ob¬ 
tained Luttinger parameters may serve as benchmark for 
theoretical models of the CU 2 O band structure. 

Conclusion. In summary, we have discovered high an¬ 
gular momentum exciton states in one-photon absorption 
spectra in high-quality CU 2 O bulk crystals. Observation 
of these states becomes possible through the rotational 
symmetry breakdown by the cubic crystal environment. 
Even though this symmetry breaking is weak due to the 
high Oh symmetry, the resulting level splittings could be 
resolved because of the for crystals record-low linewidths 
of the absorption features. 
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SUPPLEMENTARY INFORMATION 

I. SYMMETRY ANALYSIS 

In the Oh point group the P-shell excitonic states 
(quantum number I = 1) transform according to the ir¬ 
reducible representation (hereinafter we use the no¬ 
tations of Ref. [1]) whose basic functions can be chosen 
as X, y and z. The representation Vff (F-shell states, 
I = 3) of the full spherical symmetry group, ICh, is re¬ 
ducible in Oh and can be decomposed into the following 
three irreducible representations: 

= T^ + r^ + r^. (3) 

Note that in Ref. [1] the compatibility table of ICh is 
given for the Ty group, where = Fi J- r 4 -I- F 5 , so the 
following replacements are needed Fi [Ty) F^ {Oh)i 
F 4 (Td) ^ F^ {Oh) and F 5 (T^) ^ F 4 {Oh) ac¬ 
cording to the compatibility of Ty and Oh- As fol¬ 
lows from Eq. (3) the 7-fold degenerate F-shell exci¬ 
tonic state splits into the singlet F^ [basic function 
(X xyz oc Y3^2i'&,T) — Y3,-2 {'(^,t) (where f) and tp are 
the angles of r in spherical coordinates), and two triplets 
F 4 [basic functions transform like a;^, y^, , particu¬ 
larly, x^-\-iy^ oc - 3 ( 1 ?, (/?)-!-y 3 Y 3 ,i(d,v?), x^-iy^ oc 

7 ^), oc y 3 o(i?, 7 >)] and F^ 
[basic functions like x{y‘^ — z^), y{z^ — x^), z{x‘^ — y^), 
where, correspondingly, x{y‘^ — z^) -I- iy(z^ — x^) oc 
v^Y 3 , 3 (i?, (fi) - 75^3 _i(i?, (p'i x(y^ - z^) - iy(z^ - x^) oc 
-V3Y3,-3 (^,t) + \/5 Y3 ,i(i^,‘P), z(x^-y^) oc Y3^2 (^,‘P) + 
Y3-2(i?,T)]- 

With allowance for the spin and spin-orbit coupling 
the topmost valence band transforms according to the 
F^ representation and the lowest conduction band trans¬ 
forms according to the F^j" representation of the Oh point 
symmetry group, see Ref. [2] for details. For the F-shell 
envelope function one has 

P 3 - X F+ X F+ = (F^- + F 4 + F 5 ) X (F+ + F+) 

= 2Fr + F^- + 2 F 3 - + 4 F 4 + 3 F 5 . (4) 

Note, that F 4 appears 4-times, one of these representa¬ 
tions arises from the F J combination of electron and hole 
Bloch functions, namely, 

|r^) = -^|r+,-i/ 2 )|F+,i/ 2 ) 

+ ^|F+,l/2)|F+,-l/2), 
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and r 4 = Fg X FJ with basic combinations x{y^ — 
y(^^-x2)|F+) and z(a;2-y2)|r+). The F+ rep¬ 
resentation corresponds to the paraexciton Bloch func¬ 
tion, the transition to these F 4 states is spin-forbidden, 
unless those states are mixed with other states of F 4 
symmetry. 

In numerical calculations (see main text for details) it 
is convenient to consider the representation of the prod¬ 
uct of the hole Bloch function and the excitonic envelope, 
in which case we have 

Pg X F^ = (F2 -I-F4 -I-F5 ) X F^ = 2 Fg -I-F7 -|- 2 Fg , ( 5 ) 
particularly, 

F2 xFy = Fg , F4 xF^ = Fy -1-Fg , Fg xF^ = Fg -|-Fg . 

It is instructive to introduce the total momentum F of 
the hole state, where neglecting cubic symmetry effects, 
states with F = 5/2 and F = 7/2 are possible for the 
F"-shell excitons. The total angular momentum is not a 
good quantum number in the Oh point group, and the 
corresponding representations are reducible [ 11 ]: 

^ 5/2 “ ^6 + Tg , ( 6 a) 

whose basic functions are (defined as in Ref. [ 6 ]): 

( 6 b) 

|l/2), -1-1/2), 

^|-3/2) + y||5/2), ( 6 c) 

and 

^7/2 = ^6 + ^7 -I- Fg , 
with basic functions: 

Tg: ^|5/2)-i|-3/2), ^|3/2)-^| 

r.-: /Il-V 2 ) + /||l/ 2 ), 

rg- : ^|3/2) + i|-5/2), yX|-7/2)-yj|l/2), 

(7d) 


V^7/2)-Vll-l/2), i|5/2) +^1-3/2). 

The numbers in the “kets” denote the projections of F^. 
Making use of Eq. (5) one obtains the optically active 
states F 4 e Fg X F)]" and F 4 G F^ x Fg . 

The selection rules for optical transitions from the va¬ 
lence band states Eqs. ( 6 a), (7) are similar to the transi¬ 
tions in GaAs-like semiconductors (T^ point symmetry) 
with the replacement 7'T{Td) <->■ Fg (O^). 

II. COMPUTATIONAL APPROACH 

To obtain the fine structure of the F-shell exciton 
states we follow the approach of Refs. [2, 4]. We start 
in the spherical approximation, which takes into account 
explicitly the mixing of the Fy and F^ valence subbands 
and afterwards we allow for cubic contributions. 


A. Spherical approximation 

Furthermore, we neglect the central-cell corrections 
to the electron-hole potential and the electron-hole ex¬ 
change interaction. The Hamiltonian for the exciton with 
center of mass momentum P = 0 has the form [cf. Eq. 
( 1 ) of the main text]: 

where p is the momentum of relative electron-hole mo¬ 
tion, I is the angular momentum one operator acting in 
the basis of orbital hole Bloch functions F^, Sh is the 
hole spin operator (sh = 1 / 2 ), the energies are mea¬ 
sured in units of the “bulk” excitonic Rydberg TZ* = 
e^mo/( 21 i^e^ 7 (), the distances are measured in units of 
the corresponding Bohr radius, oq = h‘^e^[/{e^mo), e 
is the static dielectric constant, A = A/Rq is the di¬ 
mensionless splitting between the F^ and F^ bands, 
7( = 71 -f mo / me , 


673 -I- 472 



and 7 i (* = 1, 2, 3) are the Luttinger parameters. 

In the spherical approximation the good quantum 
numbers are F and F^, where F = L + J is the to¬ 
tal angular momentum of the envelope, L, and the hole, 
J = I + Sh, states. The wavefunction is written as 

= '^9lj{p)\LJF) 

LJ 

= E 9Lj{p)YLLAd,v)Cn:-,jjM,J^), (9) 

where gLj{p) are the radial functions, Yll^ are the spher¬ 
ical harmonics, are the Clebsch-Gordan coeffi¬ 

cients, and j J, Jz) are the basic functions of the valence 


(7a) 

-5/2), (7b) 
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band top. The effective Hamiltonian acting on the radial dices of Refs. [2, 4]. Its matrix elements between the 
functions can be constructed making use of the Appen- \LJF) and \L’J'F') states read 


:r (p" 2 ^ ... ^ 2 d, L(L + 1 ) 2 

h2- J P) = SlL’Ojj'OfF’ - ^ -^- 

\n^ r J \ ar^ r dr r 

{LJF\ . j(2)^ IL'J'F') = + 1)(2-^' + l)x 




J' F J\\J \ J‘ 

where ^ g j r the 6j-symbols defined in Ref. [5], and 


-Pl,l, 


(L-2||P<«||t) = -3S=y^i|^, 

(r + 2 ||P«||r) = P«^/ra|±5p.«,, 


L — ‘2,L 5 


with 


(f 2d L{L + 1) 


dr"^ r dr 


d^ 2L+1 d -1 


n 5 Ff — 2, L T 9 ' T ' 9 ’ 


_ d^ 2L + 1 d L(L + 2) 

PIj-\-2,L/ t 2 ^ 2 ' 

dr^ r dr r^ 


(lOa) 

(10b) 


(11a) 

(lib) 

(11c) 

( 12 ) 


In the spherical approximation the F-exciton states Making use of Eqs. (8) and (10) we arrive at the following 
correspond to total momentum F = 5/2 or F = 7/2. equations 


F = 5/2; vl; = g,^y^{p)\3, 1/2,5/2) + 5i,3/2(p)|l, 3/2, 5/2) + g3.3/2(p)|3, 3/2, 5/2), 


/ ^3.3 + y \/id-P3,l “^^-^3.3 

y^//Pi,3 (1 + m/5)-Pi.i + ^ - a -^mA.3 

V-^^^3,3 (1 - m/ 5)P3.3 + f 


\ 

A/ 


( d3.1/2\ /d3,l/2\ 

dl,3/2 1 “ I dl.3/2 j 

53,3/2/ \53,3/2/ 


P = 7/2; vi; = 53 ,i/2(p)| 3, 1/2,7/2) + 53.3/2(d)|3,3/2,7/2) + g5.3/2(d)|5,3/2,7/2), 


1 P3,3 + y 

J3M-P3.3 

-^d^^3.5 


/ 53,1/2\ 

/53,1/2\ 

^^-^3,3 

(1 — 2/r/3)P3,3 + ^ 

-A -#d^3.5 


1 53,3/2 

= A j 53,3/2 j 

\~\/Ip^5,3 

-fpP5,3 

(1 + 2/r/3)P5,5 + 1 

-aJ 

\55,3/2/ 

\55,3/2/ 


(13a) 


(13b) 


(14a) 


(14b) 


The calculated binding energies of the P-excitons and splittings between P-and P-excitons are presented in 























Fig. 3(a) and (b) of the main text. Our calculations 
show (see below for more details) that the splitting be¬ 
tween the F = 5/2 and F = 712 F-exciton states is much 
smaller compared with the splitting between the F- and 
P-excitons. Therefore the fine structure of F excitons is 
not shown in Fig. 3 of the main text. 


B. Effects of cubic symmetry 


elements of Eq. (15) can be calculated in a way similar 
to the derivation of Eq. (10b) and Eq. (A2) in Ref. [6]: 


(4) 


{L'J'F'F'^l [p(2) X ./(2) ' ' \LJFF^) 

L J m 

= 9^5 X {-l f-P'^+J+3/^x 

\/(2J + 1)(2P -b 1)(2P -b 1)(2P' -b l)x (16) 


We follow Refs. [6, 7] and include cubic contributions 
to the hole Hamiltonian in the form [cf. Eq. (2) of the 
main text] 

= 4 f E X X . 

\fe=±4 / 

(15) 

where 6 = (73 — 72 )/ 7 i- Generally, this contribution 
mixes states with different F^ for given P and with dif¬ 
ferent F. This is because in cubic symmetry the angular 
momentum is not a good quantum number. The matrix 


F' 4 P\ 
-F' m F^J 





with the 3j symbols 


a b c\ 
d e f) 


and 9j symbols 



defined in Ref. [5]. 

The wavefunctions of the states satisfy the following 
equations [12] 


P = 5/2,Fg ; 4 / = ff3,i/2(p)|3,1/2,5/2,Fg") + gi. 3 / 2 (p)|l, 3/2,5/2,Fg ) + ff3.3/2(p)|3,3/2,5/2,Fg ), (17a) 


((- 


^3,3 + f 

6 85 /e 

5 ^ 35 V 5 




(/§+-i/i) 

)Pi,3 (l + f+ ^f)P,,, + 8-A (-7S,i+?#ilP, 

+ (-¥+ + ¥' 5 ) ft.. (1 - f - fS)p,,, + 1 -a; 


1,3 


W3.1/2 

51,3/2 

V53.3/2y 


/ 53 .I /2 
= E j 5q3/2 
\53,3/2y 


(17b) 


P = 5/2, Fg ; 4- = ff3,i/2(p)|3,1/2, 5/2,Fg") + 31 . 3 / 2 ( 5 ) 11 , 3/2, 5/2, Fg ) + 33 . 3 / 2 ( 5 ) 13 , 3/2, 5/2, Fg ), (18a) 


^3,3 + y 






(1 + f-^)^^ 


40 / t) \ r> 

K I 

2 -A 


( + 5 + aZ/s) ft(,3 

(-¥+ 3,3 


1,1 


^3.1 (1 5 + 175 )^ 3.3 


Pi,3 

A 


'53.1/2 

51.3/2 

\93,3/2j 


( 53,1/2 
51,3/2 
93,3/2, 


(18b) 


P = 7/2, F/; 4 / = 33 , 1 / 2 ( 5 ) 13 ,1/2, 7/2, Ff) + 33.3/2(5)|3, 3/2, 7/2, F/) + 35.3/2(5)|5, 3/2, 7/2, F/), 


(19a) 


+3,3 + y 


(vi ^ 15 I/ 3 ) ft’.. (4^ 


3+ + 33+15) ft.3 ^ 

(+! + T03)ft.3 (l-¥-S)ft.3 + ¥A 


V- ifi 


^ + 33^ j 


_ ^+1^ _ 285 ^ p_ ^ ll _i_ _ 565 1 o , 2 


jP5.3 (1 + f -i|)P5.5 + i-Ay 


53,1/2 

53,3/2 

y95,3/2, 


' 53 . 1 / 2 ' 

= ( 53.3/2 I ■ 

\95,3/2/ 


(19b) 


(20a) 


P = 7/2, Fg ; 4 / = 33 , 1 / 2 ( 5 ) 13 ,1/2, 7/2,Fg-) + 33.3/2(5)|3, 3/2, 7/2, Fg ) + 35.3/2(5)|5, 3/2, 7/2, Fg ) 
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^ -fa,3 + f 



M _ 16\/3(5 A 

V3 35 J 

3^ 77 



fi 16V3^^ D 

(1 - ¥ + i)^’»,3 + “ - i 

- {ft‘ + ITTs) 


“ ('f»*+ ITTs) ^»'5 

(l + ¥ + i)n,6 + |-Aj 


/ ff3,l/2\ 

I 93,3/2 1 

\55,3/2/ 


/ 53,1/2\ 

E I 93,3/2 1 • (20b) 

\95,3/2 / 


F = 7/2, r-; = g3,i/2(p)|3,l/2,7/2,r8-) +g3,3/2(p)|3,3/2,7/2,r-) +g5,3/2(p)|5,3/2,7/2,r-), 


(21a) 


3,3 


16S \ 


( _ 

V V3 ^ losvsy 

(/!;■' “ 


V-+ 


231^15 


3,3 

) - ( 


(i-f- 


,* + 

2^ 


166 \ 


105V3 J 


3,3 

-2 -A 


(\/T/^+ 231^) ^3,5 ^ 


- f.- 


4.5 


3 ^ 


^)^5,3 (1+f -^)P5,5 + f-A; 


99 V5 / ^3,5 



/ 53,1/2\ 
P I 93,3/2 j 
\95,3/2/ 


I21b) 


Our calculations show that, in fact, F = 5j2 and 
F = 712 are very close to each other in energy. There¬ 
fore, we need to take into account their mixing by the 
Hamiltonian (15). Using the approach of Ref. [6] and 
Eq. (16) we obtain the following effective Hamiltonians: 




/P(5/2,r6-) dfe \ 

V 4 P(7/2,r6);’ 


(22a) 


of the four F-exciton states become almost degenerate 
and, moreover, the values of splittings between the non¬ 
degenerate states are close to the experimental data (see 
also Fig. 3(c) of the main text). The absolute values 
of the calculated transition energies deviate slightly from 
the energies measured experimentally. This deviation can 
be corrected by minor variation of the band gap energy 
in our calculations. 


where 77(5/2, Tg ) and 77(7/2, Tg ) are given by Eqs. (17b) 
and (20b), respectively, and 


+ lfPc 


= 


64 2 n 

7 \ -‘I 


^^3,3 

0 


4%/55 D \ 

■^3rf3.5\ 

0 


V+T 


3 D 4'\/55 D 

5^3,3-^P3,3 


-Mr 

33 ■ 


3,5 / 


(22b) 


for the Tg representation and 


77(rs ) ^ 


^77(5/2, Tg 




-Us 

77(7/2, Tg) 


(23a) 


where 77(5/2, Tg ) and 77(7/2, Tg ) are given by Eqs. (18b) 
and (21b), respectively, and 



/ 

0 

1 85 D 

+ 7nU5-^3,3 

45 D 

77\/4^3,5 

-^8 = 


6^/66 D 

36 -^1.3 
, 124 p 
+ ^^3,3 

0 

0 


1 

45 p 

7U3^3,3 

+ 11UT5^3,5 


(23b) 


for the Tg representation. 


III. TRANSITION ENERGIES AND 
OSCILLATOR STRENGTHS 



Figure 4: Calculated energies of the F-excitons for the princi¬ 
ple quantum number n = 4. Dots are the results of numerical 
calculations, solid lines are the description in the framework 
of the analytical model, Eq. (24). States are labeled accord¬ 
ing to the representations of the Oh group for the hole state 
Eq. (5), numbers in brackets denote the total momentum F 
of states for <5 = 0. Note that the F^ state is dark. The 
parameters of the calculation are: 77* = 87 meV, n = 0.47, 
A = 134 meV [2] 


In order to obtain more insight into the level behavior 
we introduce a two-level approximation valid for pairs of 
Tg and T^ levels for a not too large cubic anisotropy 
parameter S as follows: 


Figure 4 shows the energies of the F-excitons calcu¬ 
lated according to the developed model as function of 
the cubic anisotropy parameter <5. The other parameters 
were taken from Ref. [2]. The dots give the numerical 
results, the states are labeled according to the represen¬ 
tation of the hole state in Eq. (5). For S = —0.1 two 


fEi+ ai5 f55 \ 

135 E 2 + 0.25,) 


(24a) 


where Fi _2 are the energies of the states at J = 0, ai _2 
describe the linear-with-J shift of the diagonal energies, 
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Figure 5: Energies of E-shell states and their relative oscil¬ 
lator strengths calculated in the framework of the two-level 
model, Eqs. (24). For the F^ states two contributions to the 
oscillator strength are taken into account, see text for details. 


P describes the state mixing. The resulting energies are 


£i,2 = 


El -|- E2 + {cti + €12)5 


± 


El — E2 + {cti — 012)^ 


(24b) 


and agree well with the numerical results, see the solid 
lines in Fig. 4. 

To analyze the oscillator strengths we note that al¬ 
ready in the spherical approximation the states Tg and 
Fg originating from the F = 5/2 level are optically ac¬ 
tive. This is because the F-shell contribution with F^ 
Bloch functions is mixed with the P-shell envelope with 
Fg Bloch function, see Eq. (13a). 

The cubic Oh point symmetry of the cuprous oxide 


crystal results in mixing of the E = 5/2 and F = 1 
states described by Eqs. (22). This makes optical tran¬ 
sitions from the F^ and Fg states allowed. This is the 
only origin of optical activity for the F/" (7/2) state shown 
by the green points and lines in Figs. 4 and 5 (note 
also that the corresponding exciton wavefunction con¬ 
tains substantial contribution from the F^ Bloch exciton 
state). For the F^ states additional contributions to the 
oscillator strength arise due to the presence of the a;^, 
, z^-like envelopes of the F/" Bloch functions which, as 
a result of e.g. quartic contributions in the free-carrier 
dispersion [8] 

vt EPy +P% 

contain P-shell contributions. This contribution en¬ 
hances the oscillator strength of the r^(7/2) state as 
compared with the Fg (7/2) state. 


IV. DETERMINATION OF 73 SIGN FROM 
MICROSCOPIC BAND STRUCTURE 
CALCULATIONS 

The analysis above, see Fig. 4 and Fig. 3 from the 
main text shows that the best agreement of our calcu¬ 
lations with experiment is achieved for <5 < 0, i.e. for 
73 < 72 (we recall that 71, 72 and 73 are the Luttinger 
parameters). In order to verify this result we present in 
Fig. 6 the topmost valence band energies calculated by 
DFT in Ref. [9] and parabolic fits to these results for the 
directions [001] (left) and [III] (right). The accuracy of 
the fits is not enough to determine the values of all Lut¬ 
tinger parameters but the sign of the difference 73 — 72 
can be reliably found. The estimation yields that 73 — 72 
is negative and (73 — ^ 2 )Hi ~ —0.3 (for our parameters 
(73 - 72 )/ 7 i = H'lhi ~ -0.15). The fact that 73 < 72 
is also consistent with Ref. [10], where the DFT calcu¬ 
lations were carried out neglecting spin-orbit coupling. 
In the latter case, however, the accuracy of extraction of 
parameters is low and a direct calculation shows that in 
Ref. [10] 73 is negative. We note that our values of 71, 
72 and 73 extracted from fitting the excitonic fine struc¬ 
ture may serve as benchmark for further improvement of 
computational methods. 
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Wavevector F- >X (1 /A) Wavevector F- >R (1 /A) 

Figure 6 : Fit of the DFT band structure along the [001] and [111] axes [9]. 
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